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- - - Abstract 
00 

O I In this article, we take the point of view that the Ds{2700) be a tetraquark 

state, which consists of a scalar diquark and a vector antidiquark, and calcu- 
late its mass with the QCD sum rules. The numerical result indicates that the 
c3 , mass of the vector charmed tetraquark state is about My = (3.75 ± 0.18)GeV 

or My = (3.71 lb 0.15)GeV from different sum rules, which is about IGeV 
larger than the experimental data. Such tetraquark component should be 
very small in the Z)c,(2700). 
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cn ■ 1 Introduction 

> , 

Recently Belle Collaboration observed a new resonance 1)5(2700) in the decay 
D^D s{27Q0) D^D^K+. The resonance has the mass My = 2708±9l JoMeV, width 
Fy = 108 ± 23+i^^ MeV, and spin-parity 1-[1J. They interpret the D,(2700) as a 
^ \ cs meson, the potential model calculations predict a radially excited (cs) state 

\ with a mass about (2710 — 2720)MeV [2J. The resonance D;j(2700) is consistent 

O ' with the particle they presented previously in the 33rd international conference on 

high energy physics (ICHEP 06), My = 2715 ± lll^^MeV, Fy = 115 ± 20l^^MeV 
^ . and spin-parity 1~ [3j. In the same analysis of the DK mass distribution, Babar 

^ \ Collaboration observed a broad structure with My = 2688 ± 4 ± 3MeV and Fy = 

112±7±36MeV, which maybe the same resonance observed by Belle Collaboration 
i4j. 

In this article, we take the point of view that the vector charmed meson 1)^(2700) 
be a tetraquark state, which consists of a scalar diquark and a vector antidiquark, 
and devote to calculate its mass with the QCD sum rules [3,12]. The 1)^(2700) lies 
above the DK threshold, the decay 1)5(2700) D^K^ can take place with the 
fall-apart mechanism and it is OZI super-allowed, which can take into account the 
large width naturally. Furthermore, whether or not there exists such a tetraquark 
configuration which can result in the state Z)s(2700) is of great importance itself, 
because it provides a new opportunity for a deeper understanding of low energy 
QCD. We explore this possibility, later experimental data can confirm or reject this 
assumption. 



^E-mail, wangzgyiti@yahoo.coni.cn. 
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The article is arranged as follows: we derive the QCD sum rules for the mass of 
the Ds(2700) in section 2; in section 3, numerical result and discussion; section 4 is 
reserved for conclusion. 

2 QCD sum rules for the mass of the i^s(2700) 

In the following, we write down the two-point correlation function 11^,^ (p) in the 
QCD sum rules, 

n^.(p) = *y"rfW^-^(o|T{j^(x)j,+(o)}|o), (1) 

J^(x) = ^^'^""^ {uf{x)C'y5Cj{x)um{x)'y5lf,Csl{x) + {u ^ d)} . (2) 

We choose the vector current J/iix) which constructed from a scalar diquark and a 
vector antidiquark to interpolate the vector meson Ds{2700). 

Here we take a digression to discuss how to choose the interpolating currents for 
the tetraquark states. We can take either qq — qq type or qq — qq type currents to 
interpolate the tetraquark states, they are related to each other via Fierz transfor- 
mation both in the Dirac spinor and color space [3 E]. In this article, we take the 
qq — qq type interpolating current. 

The diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vec- 
tor C7^75, axial vector 6*7^ and tensor Ca^^. From those diquarks, we can construct 
six independent currents to interpolating the charmed tetraquark states with 1~, 



Jl{x) = ^'^""^ [u^ {.x)C-i^Cj{x)um{x)-i^^^Csl{x) + {u^ d)] , 

Jl{x) = ^^'^""^ {uJ{x)C-ff,-f5Cj{x)u^{x)^5Csl{x) + {u^d)} , 

Jf^ix) = ^^'^""^ {uJ{x)Ccj{x)u„,{x)-ff,Csl{x) + {u-^d)} , 

Jiix) = ^fe'J^'"" {uJ{x)C-ff,Cj{x)u^{x)Csl{x) + {u^d)} , 

J^ix) = ^^'^""^ {uJ{x)Caf,„Cj{x)u^{x)-f5-fuCsl{x) + {u ^ d)} , 

= ^^'J^'"" {uJ{x)C'yul5Cj{x)u„,{x)a^uCsl{x) + {u ^ d)} , (3) 



and the general current J^{x) can be written as their linear superposition, 

6 

j,{x) = (4) 

1=1 

where the Cj are some coefficients. 
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The six interpolating currents can be sorted into three types, the currents Jj^{x) 
and Jf^ix) are C75 — 6*7^75 type, the currents J^{x) and J^{x) are C — C'~^^ type, 
the currents J^{x) and are Ca^^ — 6*7,^75 type. We expect the three type 

interpolating currents result in three type of masses for the tetraquark states. 

The study with the random instanton liquid model indicates that the diquarks 
have masses about rug = 420 ± 30MeV, ttia = my = 940 ± 20MeV, = 570 ± 
20MeV PI, we expect the currents J^{x) and J^{x) interpolate the tetraquark states 
with masses larger than the ones for the currents Jj^{x) and Jf^{x). Instanton induced 
force results in strong attraction in the scalar diquark channels and strong repulsion 
in the pseudoscalar diquark channels, if the instantons manifest themselves, the 
pseudoscalar diquarks will have much larger masses than the corresponding scalar 
diquarks [lOj, the coupled Schwinger-Dyson equation and Bethe-Salpeter equation 
also indicate this fact \^\- Furthermore, the one-gluon exchange force leads to 
significant attraction between the quarks in the 0"*" channels [10]. Although the 
interpolating currents are not unique, the currents and Jf^{x) are much better 

and interpolate tetraquark states with smaller mass, we can choose either one of 
them. 

In the conventional QCD sum rules [S], there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter 
and threshold parameter sq. For the tetraquark states, if the perturbative 
terms have the main contribution, we can approximate the spectral density with the 
perturbative term. 



00 /"OO 



EmH-A / sVi7^rfs = AM^" / (5) 







where the A are some numerical coefficients, then we take the pole dominance con- 
dition, 

> 50% , (6) 



and obtain the relation. 



The superpositions of different interpolating currents can only change the contribu- 
tions from different terms in the operator product expansion, and improve conver- 
gence, they cannot change the leading behavior of the spectral density p{s) oc of 
the perturbative term. 

For the nonet light scalar mesons below IGeV, if their dominant Fock compo- 
nents are tetraquark states, even we choose special superposition of different currents 
to weaken the contributions from the vacuum condensates to warrant the main con- 
tribution from the perturbative term, we cannot choose very small Borel parameter 
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to enhance the pole term. For small enough Borel parameter M^, the per- 
turbative corrections of order 0{as{M)), C(a^(M)), maybe large enough to 
invalidate the operator product expansion. 

We can choose the typical energy scale /i = M = IGeV, in that energy scale, 
y/sQ ~ 2.2GeV. There are many scalar mesons below 2GeV [12], their contributions 
are already included in at the phenomenological side. Pole dominance cannot be 
fully satisfied for the tetraquark states with light flavor. 

Failure of pole dominance do not mean non-existence of the tetraquark states, 
it just means that the QCD sum rules, as one of the QCD models, may have short- 
comings. We release some criteria and take more phenomenological analysis, i.e. we 
choose larger Borel parameter to warrant convergence of the operator product 
expansion and take a phenomenological cut off to avoid possible comminations from 
the high resonances and continuum states |13j . 

If we insist on to retain pole dominance besides convergence of the operator 
product expansion in the QCD sum rules for the tetraquark states, the hidden 
charmed and bottomed tetraquark states, and open bottomed tetraquark states 
may satisfy the criterion in Eq.(7), as they always have larger Borel parameter 
and threshold parameter sq. 

For examples, in Ref.[Tl], the authors take the X(3872) as hidden charmed 
tetraquark state and calculate its mass with the QCD sum rules, the Borel parameter 
and threshold parameter are taken as = (2.0 — 2.8)GeV^ and sq = (17— 18)GeV^; 
in Ref. [15j. the authors take the Z(4430) as hidden charmed tetraquark state and 
calculate its mass with M2 = (2.5-3.1)GeV2 and Sq = (23-25)GeV^ furthermore, 
the authors calculate the corresponding bottomed one, and choose = (8.0 — 



8.3)GeV2 (or = (8.0 - 9.9)GeV2) and sq = 125GeV^(or sq = 135GeV2). In 



those sum rules, although the windows for the Borel parameters are rather small, the 
as{M) is small enough to warrant convergence of the operator product expansion, 
the relation in Eq.(7) can be well satisfied. 

The correlation function Il^^{p) can be decomposed as 



due to Lorentz covariance. The invariant functions Hi and IIq stand for the contri- 
butions from the vector and scalar mesons, respectively. In this article, we choose 
the tensor structure g^u — to study the mass of the vector meson. 

According to basic assumption of current-hadron duality in the QCD sum rules 
[5], we insert a complete series of intermediate states satisfying unitarity principle 
with the same quantum numbers as the current operator J^{x) into the correlation 
function Il^yip) to obtain the hadronic representation. After isolating the pole term 
of the lowest state Ds(2700), we obtain the following result: 
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where we have used the following definition, 



(0|J^(0)|D,(2700)) = fvM^e^, (10) 

here is the polarization vector of the Z^s(2700) and fy is the residue of the pole. 

In the following, we briefiy outline operator product expansion for the correlation 
function n^j,(p) in perturbative QCD theory. The calculations are performed at 
large space- like momentum region <^ 0, which corresponds to small distance 
a; ~ required by validity of operator product expansion. We write down the "full" 
propagators Sij{x) (the Uij{x) and Dij{x) for the u and d quarks can be obtained with 
a simple replacement of the nonperturbative parameters) and Cij{x) of a massive 
quark in the presence of the vacuum condensates firstly jSj^ . 

192 ^ ^ ' 1152 ^ ^ 

-^G%{^a^'' + a^-" ^) + ■ ■ ■ , (11) 



32'k'^x 



(27r)4 J \f^-m^ 4 (k'^ - m. 



asGG P + m^U , 

+y(^)^.-^(p3^ + ---^ (12) 



where {sQsCtGs) = {sggCrafjG'^^ s) and (^a^) = ( ""'^"^'^ — ), then contract the quark 
fields in the correlation function n^j,(p) with Wick theorem, and obtain the result: 

n^i^(p) = i^kijek'i'j'ekmnek'm'n' j d^Xe'^'^'Tr {'j5'Jf,CS^,„{-x)G'Ju75Um'mi-x)} 

Tr {^,G,Axh,CU^A^)C} . (13) 

Substitute the full s, c and u quark propagators into above correlation function and 
complete the integral in coordinate space, then integrate over the variable k, we can 
obtain the correlation function ni(p^) at the level of quark-gluon degrees of freedom: 



^One can consult the last article of Ref. 5 for technical details in deriving the full propagator. 
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log/C 



rUsP 
192^ 
mdqq) 



[ dt [6{t - l){qq) + {t'^ + t-2) (ss)] /C log /C 

^0 



1927r4 
rus 



1287r4 
mc{qgs(rGq) 



1287r4 



+ i)/C log/C 



log/C 



3847r4 

^2 



TUsP 

384^ 

{qq){ss) 



mcrris 
' 247r2 



[(3t + 1) (s^sf^Gs) + 12(g53CTG'g)] /C log /C 
c^t [(t - t^){sgsaGs) + 6(t - t^)(g^5(TG'g)] log/C 
{m){ss)p' ^'dt(t-t^) log/C 



di/Clog/C+ ^27^2 
dt[2(^?)2 + i(^g)(ss)] log/C 



{qq){sgs(jGs) + {ss){qgsaGq) 
247r2 



^0 



dtt log/C , 



(14) 



where /C(p2) = (1 - t)m2 - t(l - t)p^. 

We carry out operator product expansion to the vacuum condensates adding up 
to dimension-8. In calculation, we take assumption of vacuum saturation for high 
dimension vacuum condensates, they are always factorized to lower condensates with 
vacuum saturation in the QCD sum rules, factorization works well in large Nc limit. 
In this article, we take into account the contributions from the quark condensates, 
mixed condensates, and neglect the contributions from the gluon condensate. In 
calculation, we observe the contributions from the gluon condensate are suppressed 
by large denominators and would not play any significant roles. 

Once analytical results are obtained, then we can take current-hadron duality 
below the threshold sq and perform Borcl transformation with respect to the variable 

— —p^i finally we obtain the following sum rules: 
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Ml 
M2 



Ml 



"0 , Imn(s) 
as exp 

"° , Imn(s) 
as sexp 

, Imn(s) 
as exp 

,2 TT 



M2 
s 



(15) 



(16) 



Imn(s) 



TT 



614407r6 



dt 



1927r4 £ [^^^ ~ ^^^^^^ + (t' + t - 2)(ss)] /C 

"^c(gg) 



19277^ 



+ 



1287r4 



TUciggscrGq) 



^{qq) + (l + ^){ss) 



dt{l + -)}C 

A t 



3847r4 



TUgS 

"384^ 



dt [{3t + l){sgsaGs) + 12{qg,aGq)] /C 
rft [(t - t^){sgsaGs) + 6(t - t^)(g^sO-G'g)] 



rricms 



dtJC - 



{qq){ss)s 



1 

dt[2{qqy + t{qq){ss)] 



[ dt{t-f 

J A 



247r2 



{qq){sgsCrGs) + {ss) (qgsCrGq ) 
247r2 



dtt . 



(17) 



where A 



3 Numerical result and discussion 

The input parameters are taken to be the standard values (qq) = — (0.24±0.01GeV)'^, 
{ss) = (0.8 ± 0.2) (gg), {qgscrGq) = ml{qq), (sgsCrGs) = mg(ss), = (0.8 ± 
0.2)GeV^ m, = (0.14 ± 0.01)GeV and = (1.4 ± 0.1)GeV il [TH]. For the 
multiquark states, the contribution from terms with the gluon condensate (^^^) is 
of minor importance [13], and the contribution from the (^^^) is neglected here. 
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perturbative term 


+96% 


(gg), {ss) 


+33% 


(qgsCxGq), {sgsCrGs) 


-10% 


{qqy,{qq){ss) 


-24% 


{qq){sgs(rGs), {qgsaGq){ss) 


+4% 



Table 1: The contributions from different terms in Eq.(15) for Sq = 16GeV^ and 
M2 = 6GeVl 

From Table 1, we can see that the dominating contribution comes from the 
perturbative term, (a piece of) standard criterion of the QCD sum rules can be 
satisfied. If we change the Borel parameter in the interval = (5 — 7)GeV^, the 
contributions from different terms change slightly. 

Although the contributions from the terms concerning the quark condensates 
and mixed condensates are rather large, however, they are canceled out with each 
other, the net contributions are of minor importance. Which is in contrast to the 
sum rules with other interpolating currents constructed from the multiquark con- 
figurations, where the contribution comes from the perturbative term is very small 
|17j . the main contributions come from the terms with the quark condensates (qq) 
and (ss), sometimes the mixed condensates (qgsO'Gq) and {sggCrGs) also play im- 
portant roles (for example, the first three articles of the Ref.|13]). One can choose 
special superposition of different currents to weaken the contributions from the vac- 
uum condensates to warrant the main contribution from the perturbative term, it 
is somewhat of fine-tuning. 

The values of the vacuum condensates have been updated with the experimental 
data for r decays, the QCD sum rules for the baryon masses and analysis of the 
charmonium spectrum [TB]. As the main contribution comes from the perturbative 
term, uncertainties of the vacuum condensates can only result in very small uncer- 
tainty for numerical value of the mass My, the standard values and updated values 
of the vacuum condensates can only lead to results of minor difference, we choose 
the standard values of the vacuum condensates in the calculation. 

In Fig.l, we plot the value of the My with variations of the threshold parameter 
So and Borel parameter M^. If ^/sq < 3.55GeV, My > sq, we cannot take into 
account all contributions from the Ds(2700), furthermore, the My changes quickly 
with the variation of the Borel parameter M^, the threshold parameter sq should be 
chosen to be ^/sq > 3.6GeV. The value of the My is almost independent on the Borel 
parameter at about ^/so = 4.0GeV. In this article, the threshold parameter sq 
is chosen to be sq = (16 ± 2)GeV^. It is large enough for the Breit-Wigner mass 
My = 2708 ± g^JjMeV, width Ty = 108 ± 23^11 MeV. However, the standard 
criterion of pole dominance cannot be satisfied, the contribution from the pole term 
is less than 13%. If one insist on that the multiquark states should satisfy the same 
criteria as the conventional mesons and baryons, the QCD sum rules for the (light 
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Figure 1: My with Borel parameter and threshold parameter sq. 

and charmed) tetraquark states should be discarded. For detailed discussions about 
how to select the Borel parameters and threshold parameters for the multiquark 
states, one can consult Ref . [13] . 

Taking into account all the uncertainties, we obtain the value of the mass of the 
Ds{2700), which is shown in Fig.2, 

My = (3.75±0.18)GeV. (18) 

It is obvious that our numerical value is larger than the experimental data My = 
2.708GeV, the vector current can interpolate a charmed tetraquark state with the 
mass about My = 3.75GeV or even larger, such tetraquark component should be 
small in the ^,(2700). 

If one want to retain the pole dominance of the conventional QCD sum rules, we 
take the replacement for the weight functions in Eqs. (15-16), 
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Figure 2: My with Borel parameter from Eq.(16). 
and obtain new QCD sum rules for the mass of the vector tetraquark state. 




As the main contributions come from the perturbative term, the hadronic spec- 
tral density above and below the threshold can be successfully approximated by the 
perturbative term. If we take typical values for the parameters ^/so = 4.0GeV and 
= (7 — 9)GeV^, the contribution from pole term in Eq.(20) is dominating, about 
53% — 84%. Taking into account all the uncertainties, we obtain the value of the 
mass of Ds(2700), which is shown in Fig. 3, 

Mv = (3.71 ±0.15)GeV. (22) 



4 Conclusion 

In this article, we take the point of view that the Ds{2700) be a tetraquark state 
which consists of a scalar diquark and a vector antidiquark, and calculate its mass 
with the QCD sum rules. The numerical result indicates that the mass of vec- 
tor charmed tetraquark state is about My = (3.75 ± 0.18)GeV or My = (3.71 ± 
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Figure 3: My with Borel parameter from Eq.(21). 

0.15)GeV, which is about IGeV larger than the experimental data. Such tetraquark 
component should be very small in the D<j(2700), the dominating component may 
be the cs state, we can take up the method developed in Ref.[TB] to study the mixing 
between the two-quark component and tetraquark component with the interpolating 
current J^{x) = cos9J^{x) + sin9{qq)s{x)'j^c{x). The decay 1)^(2700) D^K^ can 
occur mainly through creation of the uu pair in the QCD vacuum, we resort to the 
^Po model to calculate the decay width [19], although the ^Pq model is rather crude. 
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